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Within the classical field model, we find that the phase of a Bose-Einstein condensate undergoes 
a true diffusive motion in the microcanonical ensemble, the variance of the condensate phase change 
between time zero and time t growing linearly in t. The phase diffusion coefficient obeys a simple 
scaling law in the double thermodynamic and Bogoliubov limit. We construct an approximate cal- 
culation of the diffusion coefficient, in fair agreement with the numerical results over the considered 
temperature range, and we extend this approximate calculation to the quantum field. 
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I. INTRODUCTION 

Phase coherence is one of the most prominent prop- 
erties of Bose-Einstein condensates, relevant for apphca- 
tions of condensates in metrology and quantum informa- 
tion The issue of condensate phase dynamics and 
phase spreading at zero temperature due to interactions 
has been extensively studied in theory Q and experi- 
ments 1^, i^] . There is a renewed interest in this issue 
of temporal phase coherence due to the recent studies 
in low dimensional quasi condensates, both experimen- 
tally d, 0, [l] and theoretically [9]. The present work 
addresses the problem of the determination of fundamen- 
tal limits of phase coherence in a true three dimensional 
Bose-Einstein condensates at non zero temperature. 

The effect of finite temperature on phase coherence in 
a Josephson junction realized by a condensate trapped 
in a double well potential has been studied in and 
[m . The situation is different when the two condensates 
are separated. In this case there is no restoring force for 
the relative phase which then evolves independently in 
the two BECs The effect of the temperature in this 
case, joint to the effect of the interactions, is to provide 
a spreading in time of the relative phase. 

In a previous work [l^, we considered a condensate 
prepared in an equilibrium state in the canonical ensem- 
ble. In that case we could show using ergodicity that the 
phase change of the condensate during a time t has a vari- 
ance which grows proportionally to . In other words, 
the condensatephase spreading in the canonical ensem- 
ble is baUistic and not diffusive [H, [11, [13, El] ■ As we 
could calculate in [ll| using an ergodic theory, the coeffi- 
cient of this super-diffusive thermal spreading is propor- 
tional to the variance of the energy in the considered equi- 
librium state. If we now suppress the fluctuations of en- 
ergy in the initial state, by moving from the canonical en- 
semble to the microcanonical ensemble, the ballistic ther- 
mal spreading disappears and one may expect that the 
condensate phase undergoes a genuine diffusion in time. 
In the present work we show that this is indeed the case 
and we study this genuine phase diffusion numerically 
within the classical field model [H [13, [U, [H, [H, H] 



described in section |TT1 The numerical results are pre- 
sented in section IIIIl and their analysis shows the exis- 
tence of simple scaling laws and of an universal curve 
giving the phase diffusion coefficient in the double ther- 
modynamical limit {N — > oo, volume V — > cxd density 
p = constant) and Bogoliubov limit (TV — > oo, coupling 
constant 5 — > 0, Ng = constant). In section [TVl we derive 
an aproximate formula for the diffusion coefficient, that 
we compare to the numerical results and that we also 
extend to the quantum field. We conclude in section IVl 



II. CLASSICAL FIELD MODEL AND 
NUMERICAL PROCEDURE 

We consider a lattice model for a classical field ■!/;(r) 
in three dimensions. The lattice spacings are l\, l2, I3 
along the three directions of space and dV = hhh is 
the volume of the unit cell in the lattice. We enclose 
the atomic field in a spatial box of sizes Li, L2, and 
volume V = L1L2L3, with periodic boundary conditions. 
To guarantee efficient ergodicity in the system we choose 
non commensurable square lengths in the ratio : L2 ■ 
L§ = \/2 : (1 + a/5)/ 2 : VS. The lattice spacings squared 
If, I2, are in the same ratio. 

The field ^ may be expanded over the plane waves 



V 



(1) 



where k is restricted to the first Brillouin zone, ka G 
[— tt/Zq, tt/Zq] and a labels the directions of space. 

We assume that, in the real physical system, the total 
number of atoms is fixed, equal to N. In the classical 
field model, this fixes the norm squared of the field: 



dF^|^(r)|2=7V. 



Equivalently the density of the system 
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is fixed for each realization of the field. The evolution of 
the field is governed by the Hamiltonian 

H = Y, EkaW + f E dVr{v)r (r)V'(r)V'(r), (4) 



where is the dispersion relation of the non-interacting 
waves, and the binary interaction between particles in the 
real gas is reflected in the classical field model by a field 
self-interaction with a coupling constant 



(5) 



where a is the s-wave scattering length of two atoms. 

As a matter of a fact we use here the same refinement 
as in consisting in modifying the dispersion relation 
in order to obtain for the ideal gas the correct quantum 
values of the mean occupation numbers at equipartition 
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However we do not expect this to have a large impact 
here as we put a cut-off at an energy of the order of fcsT. 
More precisely we choose the number of the lattice points 
in a temperature dependent way, such that the maximal 
Bogoliubov energy on the lattice is equal to /cbT: 

maxk [{h^k^ /2m){2pg + h^k^ /2m)YI'^ = keT . (7) 

The discretized field has the following Poisson brackets 



z;i{V'(ri),^*(r2)}- 



dV 



(8) 



where the Poisson brackets are such that df /dt — {f, H} 
for a time-independent functional / of the field ^. The 
field then evolves according to the non linear equation 

M 



parametrize the microcanonical ensemble by the temper- 
ature T such that the mean energy of the field in the 
canonical ensemble at temperature T is equal to E. 

To generate the stochastic initial values of the clas- 
sical field we proceed as follows, (i) First we generate 
1000 stochastic fields in the canonical ensemble at tem- 
perature T, as explained in [13], and we compute the 
average energy of the field (-E)can and its root mean 
squared fluctuations a — VVar E. (ii) We generate other 
fields, still in the canonical ensemble, and we filter them 
keeping only realizations with an energy E such that 
\E - (-B)canl < 0.01 a/2, (ni) We let each field evolve for 
some time interval with the Eq. ^ to eliminate transients 
due to the fact that the Bogoliubov approximation, used 
in the sampling, does not produce an exactly stationary 
distribution. After this 'thermalization' period we start 
calculating the relevant observables, as ip evolves with 
the same equation ([9]). In practice this equation is in- 
tegrated numerically with the FFT splitting technique. 
The ensemble of data reported here has required a CPU 
time of about two years on Intel Xeon Quad Core 3 GHz 
processors. 



III. NUMERICAL RESULTS AND SCALING 
LAWS 

The first important result that we obtain is the diffu- 
sive behavior of the condensate phase. In Fig[T]we show 
an example of numerical data. From bottom to top, five 
values of ksT/pg are presented for a constant number of 
atoms N — 2.36 x 10^. The wavy line with error bars is 
the phase variance as a function of time obtained with 
about 1200 stochastic realizations 0]. The solid line is a 
linear fit from which we deduce the value of the diffusion 
coefficient. 
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ih dtip — < ksT 
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+ 5|V(r,t)P[> V- 
(9) 



We introduce the density and the phase of the conden- 
sate mode 



ao 



(10) 



The quantity of interest is the variance of the condensate 
phase change during t: 



where 



^{t) = e{t) - 0{Q). 



(11) 



(12) 



The averages are taken over stochastic realizations of the 
classical field, as the initial field samples the microcanon- 
ical ensemble with an energy E. For convenience, we 



The diffusive behavior of the condensate phase is 
strictly related to the long time behavior of the time cor- 
relation function C of the condensate phase derivative ip, 

c{\t' - = m'm")) ~ m'))m")) , (w) 

where we used the fact that C depends only on \t' — t"\ 
for a steady state classical field. By writing (p{t) in terms 
of its time derivative, one obtains [131 



Var(p(t) = 2t dr C(t) - 2 dr tC{t) 



(15) 



If C{t) has a non-zero limit at long times, as it was the 
case in the canonical ensemble [l^ , Var f grows quadrati- 
cally in time. Here, in the microcanonical ensemble Var ip 
grows linearly in time and we expect that C{t) when 
i — > oo. An illustration of that, for two values of the 
temperature, is given in Fig[2] where, for convenience. 
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FIG. 1: Variance of the condensate phase change tfi(t) as a 
function of time. Wavy hne with error bars: Numerical re- 
sults. Solid lines: A linear fit. From bottom to top, the 
reduced temperature ksT/pg is 9.7, 13.2, 15.7, 19.6, 24.2. 
The number of atoms is fixed to A*' = 2.37 x lO". The 
high energy cut-off is fixed according to ([7}, on a grid 32^, 
so that the temperature slightly varies, from bottom to top: 
kBT/in^ /mV'^'^) = 16864,16411,16212,16010,15854. The 
time is in units of mV^^^ /h. 



C{t) is calculated with a simplified formula for the phase 
derivative |27|| 



Tup ~ -pg - ^^{Uk 



14)' l&K 



(16) 



In p6)) the 6k; are the field amplitudes on the Bogoliubov 
modes 

We now investigate numerically how the diffusion co- 
efficient scales in different limits. First we consider the 
"Bogoliubov limit" introduced in [1^ 



N oo , g ^ with Ng = constant , 



(17) 



the other parameters (V^, li, ^2, ^3, ?") being fixed [30|. In 
this limit, the number of non condensed particles con- 
verges to a non zero value while the non condensed frac- 
tion vanishes. The time evolution of the Bogoliubov 
occupation numbers rik = |6kP is then mainly due to 
terms in the interaction Hamiltonian which are cubic in 
the non condensed field amplitude and linear in the con- 
densate amplitude and thus of order e = gVN. Physi- 
cally these cubic terms describe interactions among Bo- 
goliubov modes such as Landau and Beliaev processes 
, which are included in the classical field 
model [la tia M m m, HE [13, m, [M lU . They lead 
to evolution rates of the nk of order e^. We thus expect a 
phase diffusion coefficient of the same order , which is 
(X 1/N according to (fT7|) . This expectation is confirmed 
numerically as we show in FiglHl where we find that DN 
is constant within the error bars over a factor 5 variation 
of and for three considered temperatures. 

We now investigate the existence of a thermodynamical 
limit for the quantity DN, given that the Bogoliubov 
limit is already reached. The thermodynamical limit is 




t[mV" /ii] 

FIG. 2: Correlation function C{t) of the phase derivative (p{t) 
given by the non-oscillating approximation (|16|) , as a function 
of time. Solid line: Numerical results. Dashed line: Result of 
Bogoliubov theory. Dashed-dotted line: Prediction of the pro- 
jected Gaussian approach of section HVl The number of atoms 
is fixed to A = 5 X lO'' , and the Gross-Pitaveskii chemical po- 
tential is fixed to pg — 700h^ /mV^^^ . In (a) the temperature 
is ksT/ih^/mV'^^^) = 5469, with a grid size 18^ In (b) the 
temperature is ksT /{h? /mV'^^^) — 14054, with a grid size 
30^ 



defined as usual as 



A^ — > oo , y — > oo with p ~ constant , 



(18) 



the other parameters {g,li,l2,l3,T) being fixed. The re- 
sult is shown in FigH] where DN is constant within the 
error bars, over a factor 5 of variation of N and 4 con- 
sidered temperatures. 

In what follows, using dimensional analysis, we show 
that for our cut-off procedure (O, the dimensionless 
quantity HDN/pg is a function of a single parameter 
ksT/ pg once the Bogoliubov and thermodynamical lim- 
its are reached. Six independent physical quantities are 
present in the model 



5, V, keT, N}. 



(19) 



The lattice spacings hjh, h are not independent param- 
eters since their ratios are fixed and their value is deter- 
mined by ([7]) once the quantities are fixed. Equiva- 
lently we can replace g by pg and the volume V by ksTc, 
where Tc is the transition temperature of the ideal gas 
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FIG. 3: Scaling of the phase diffusion coefficient in the Bo- 
goliubov limit (|17|l . for a factor 5 variation of the atom num- 
ber A'^. The Gross-Pitaevskii chemical potential is fixed to 
pg = 700h^ /mV^^^ . Points with error bars: Simulation re- 
sults. The lines connect the points with the same temper- 
ature. From bottom to top: ksT = 5469h^ /mV'^^^ with a 
grid size 18^, ksT = 6606h'^ /mV^^^ with a grid size 20^, 
ksT = 9231h'^ /mV^^'-^ with a grid size 24^ 
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FIG. 4: Scaling of the phase diffusion coefficient in the ther- 
modynamic limit (|18p . for a factor 5 variation of the atom 
number A'^. Points with error bars: Simulation results. The 
lines connect the points with the same kBT/pg. The points 
most on the left of the figure, with A'^ = 2.37x 10® , are the ones 
of Fig[T] with a grid size 32"^ . The other points are for a grid 
size 24^ {N = 10^) and for a grid size 18^ [N = 4.22 x 10^). 



given by 



We then have 
hDN 



rnksTc 



3/2 



- C(3/2) 



pg 



pg pg . 



(20) 



(21) 



The three quantities h, m, pg can be recombined to form 
a length, a time and a mass which are three indepen- 
dent dimensioned quantities. Since / and its other three 
variables are dimensionless, / does not depend of its first 
three variables. In the thermodynamical limit iV — > oo 



so the sixth variable of / drops out of the problem. In 
the Bogoliubov limit, kBTd pg ^ oo so that the forth 
variable of / also drops. We thus conclude that 



hDN _ ^ fkgT 



pg 



pg 



(22) 



In FigOwe show the graph of / as obtained by our clas- 
sical field model collecting all the simulation results of 
Fig|4]and FigO We have used a log- log scale in Fig[5]to 
reveal that the function / is approximately a power law 
in the considered range of ksT/pg. 
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FIG. 5: Universal curve for the rescaled phase diffusion co- 
efficient as a function of ksT/pg in log-log scale. Symbols 
with error bars: Simulation results. Circles: Results of FigO 
Squares: Results of FigU The solid line connects the points 
with the largest value of A'' of Fig|4] plus the average of the 
two points of FigO with the lowest temperature. Dashed- 
dotted line with filled diamonds: D^pp'^°^ from the projected 
Gaussian approach. The Fk for the projected Gaussian ap- 
proach are calculated on the same discrete grids used in the 
simulation points connected by the solid line. 



IV. PROJECTED GAUSSIAN 
APPROXIMATION 

In this section we propose an approximate analytical 
formula for the phase diffusion coefficient that gives some 
physical insight and can be extended to the quantum field 
case. 



A. Classical field 

We wish to calculate the integral of a correlation func- 
tion C{t) = {A{t)A{0)) - (A(i)) (yl(0)) of an observable A 
of the form 



(23) 



where 6k are the amplitudes of the field over the Bogoli- 
ubov modes. We thus introduce the M. M. covariance 
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matrix Q with matrix elements 

Qk,k'(i) = ('5nk(t)<5nk'(0)) 



(24) 



where 5n\^ = rik — nvi is the fluctuation of the occupation 
number of the corresponding Bogoliubov mode and Ai = 
V/{lxl2h) — 1 is the number of Bogoliubov modes. One 
thus has 



C{t)= A-Q{t)A 



(25) 



where A is the vector of components Ak. Since the sys- 
tem is described in this section by the microcanonical 
ensemble for the Bogoliubov Hamiltonian, the matrix Q 
obeys the relation 



?-Q = *0 and g?=0, 



(26) 



where the vector e* collects the Bogoliubov energies. 

By using the microcanonical classical averages [i^ one 
directly accesses the t = value of the matrix Q, 



for k 7^ k' 
for k = k' 



Qk,k'(0) 



M + 1 



gk,k(0) - ^ 



(27) 
(28) 



Remarkably we can express this result in terms of the 
result one would have in the canonical ensemble with av- 
erage energy equal to the microcanonical energy, adding 
a projector which suppresses energy fluctuations: 



Q(0) 



M 



M 



ptgGauss^Q^p 



with 



Pi 



k.k' 



k,k' 



Ck^k' • 



(29) 



(30) 



The vector a is adjoint to the vector e* so that Pe ~ Q. 
Its components are given by 



1 



ak 



Mek 



(31) 



and QGauss value of the covariance matrix in the 

canonical ensemble 



")Gauss 
Sk,k' 



(0) 



3k,k"^k ■ 



(32) 



The apex "Gauss" reminds the fact that the 6k have a 
Gaussian probability distribution in the canonical ensem- 
ble, contrarily to the case of the microcanonical ensemble. 

For the t = value of the phase derivative correlation 
function one then obtains 



C(0) 



with 
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M + 1 



kj^O 



k"fe 
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A-^Uk 



A^ = -^{Uk + Vkf ■ 



(33) 



(34) 



We verified that ([33]), represented as a dashed line in 
FigEl is in agreement with the numerical simulation as 
one enters the Bogoliubov limit (fT7|) . Note that within 
the Bogoliubov aproximation b\^{t) ~ 6k(0) e""''*'', the 
phase derivative correlation function remains equal to its 
< = value (l33|) at all times. This is in clear disagreement 
with the numerical simulation, and it would lead to a 
ballistic spreading of the condensate phase. 

Our approximate treatment consists in extending the 
relation ((29|l at positive times, using the fact that in a 
Gaussian theory one would have 



Q^r^(i) = Qk,k'''(0)e 



Gauss / 



-,-rkt 



One indeed assumes in the Gaussian model 



l(&k(<)^'k'(0)) 



Gauss 



k.k'^ifce 



-Tkt/2 



(35) 



(36) 



and one uses Wick theorem to obtain ([55)1 . Physically 
equation ((36)) describes Landau Beliaev processes that 
decorrelate the 5k. It can be derived for example with a 
master equation approach as done in an appendix of [T^] . 

For the phase derivative correlation function one then 
obtains the approximate expression 



^approX(-^>) 
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k'^fc'e ^i-'* 



(37) 



We represent ([37| as a dashed-dotted line in FiglS] The 
resulting approximation on the diffusion coefficient is ob- 
tained by integration 



p+oo 

pjapprox^ / C''PP'-™(t)dt. 

Jo 



(38) 



In FigO we compare the approximation (|38p (diamonds 
linked by a dashed-dotted line) to the numerical simu- 
lation results. The agreement is acceptable in the con- 
sidered range of ksT/ pg. The Landau-Beliaev damping 
rates Fk are calculated on the same discrete grid as the 
simulation points as explained in [4l|. 



B. Quantum field 

In this subsection we extend the approximate formula 
for the phase diffusion coefficient to the quantum case. 
The Bogoliubov amplitudes and occupation numbers are 
now operators 6k, "-k- As in the classical case we intro- 
duce the covariance matrix of the Bogoliubov occupation 
numbers 



Qk,k'(t) = (<5nk(i)'5rik'(0)). 



(39) 
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To obtain the t = value of Q we need to compute 
quantum averages in the microcanonical ensemble. To 
this end we use a result derived in [l^] giving the first 
deviation between the microcanonical expectation value 
(O) and the canonical one (O)can(r) in the limit of a 
large system for an arbitrary observable O: 

(40) 

where the canonical temperature T is such that the mean 
energy (-ff)can in the canonical ensemble is equal to the 
microcanonical energy E. 

For k ^ k', using dhk/dT — ekfikink + l)/fcsr^, one 
finds 



(9k,k'(0) 



f-k(-k'nk{nk + l)nk'{nk' + f) 

Eq^Oeg%(?\ + 1) 



(41) 



This scales as in the thermodynamic limit. Since 

the number of off-diagonal terms of Q in (P5)l is about 
M times larger than the number of diagonal terms of Q, 
we have for consistency to calculate the diagnal terms 
of Q up to order that is the deviation from the 

canonical value nk{nk + 1) is not required. To exactly 
obtain the energy conservation (j26p . it is however conve- 
nient to include, rather than the exact deviation between 
the canonical and microcanonical values, an ad hoc ap- 
proximate correction of order 1 /A4 : 



QkM ^ nkiuk + 1) - 



Eq^O + 1) 



(42) 



In this way, we recover the structure of where the 
t = value of Q is deduced from the one in the canonical 
ensemble. 



Q^x''{t^O)^nk{nk + l)Si,M' 
by the action of a projector P, 

Q{t === 0) ~ P^Q^''''''{t = 0)P. 



(43) 



(44) 



The projector P still involves the dyadic structure ([30)) . 
with a new expression for the vector a: 



eknkjnk + 1) 



(45) 



As a check, one can apply the classical field limit to the 
above quantum expressions. One recovers (|29p . apart 
from the global factor A4/{A4 -1-1), whose deviation from 
unity gives rise to terms beyond the accuracy of the 
present calculation. 

At positive times, our quantum projected Gaussian ap- 
proximation assumes that (|44)) still holds. 



(46) 



gapprox^^^ _ ptgGauss^^^p 

with the Gaussian covariance matrix 

Qkr'(i) = Sy^M'Tikiuk + l)e 



(47) 



where the Landau-Beliaev damping rate Fk is now the 
usual one, that is for the quantum field theory. From 
()25|) one obtains an approximate expression for the phase 
derivative correlation function, and from (|38p an approx- 
imate expression for the quantum field phase diffusion 
coefficient [4211: 



^approx ^ J2[{PA] 



• nkjiik + 1) 



(48) 



where the projection of the vector A was introduced: 



(Uk + Vkf - c. ggl^g + + 1) 



(49) 



In this expression, the modes amplitudes Uk , Vk and en- 
ergy ek have the usual expressions of the quantum field 
Bogoliubov theory. 



Uk + Vk = 



1 



Uk-Vk 



2pg + h^k'^/2m 



1/4 



2m 



2W 



2m 



1/2 



(50) 
(51) 



In appendix [A] we give the explicit expression of the 
Landau-Beliaev damping rates F^ and the approximate 
phase diffusion coefficient in the thermodynamical limit. 
The existence of such a limit is due to the fact that none 
of the momentum integrals involved are infrared diver- 
gent, keeping in mind that {Uk + Vk)'^, ek and F^ vanish 
linearly with k [s^, while nk{nk + 1) diverges as 1/k^. 
As expected from the analysis of the previous section, the 
scaled diffusion coefficient fiD'^PP'^°^N/ pg is a function of 
ksT/pg only, see (|A16I) . 

In Figiniwe show the values of /J^pprox ^j^g thermo- 
dynamic limit for the same values of ksT/ pg as in Fig|5l 
To see the effect of an energy truncation at /csT, we also 
show the values of /japprox obtained by introducing an 
energy cut-off €q < fc^T in (|A16p and the same cut-off 
Efc < fc^r in the integrals (jAllAOp giving the damping 
rate Fg. As expected, the resulting values of £)approx ^^.^ 
close to the values of D^-pprox obtained for the classical 
field model in FigO and that we have reported in FiglH] 
for comparison. We conclude that the diffusion coefficient 
is indeed affected by an energy cut-off, and the coefficient 
obtained in the classical field simulations with an energy 
cut-off ksT might differ quantitatively from the real one 
by a factor of about two. 



V. CONCLUSION 

Using a classical field model, we have shown that the 
phase of a Bose-Einstein condensate undergoes true diffu- 
sion in time, when the gas is initially prepared in the mi- 
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FIG. 6: Approximate phase diffusion coefRcient in the ther- 
modynamic hmit as a function of ksT/pg in log-log scale. 
Solid line with crosses: d^pp'^°^ for the quantum field from 
Ea. (IA16p . Dotted line with crosses: Effect of an energy cut- 
off equal to ksT in Ea. HA16|l . Dashed-dotted line with filled 
diamonds: D=»pp'^°^ for the classical field from Ea. (f38|l (same 
data as in Fig[5]). 



crocanonical ensemble. Parametrizing the microcanoni- 
cal energy E by the temperature T of the canonical en- 
semble with average energy E, we could show that the 
rescaled diffusion coefRcient hDN/pg, where N is the 
fixed number of particles and pg is the Gross-Pitaevskii 
chemical potential, is a function of a single variable 
ksT/pg in the double thermodynamic and Bogoliubov 
limit. 

We have derived an approximate formula for the diffu- 
sion coefficient, in fair agreement with the classical field 
simulations. We could generalize the approximate for- 
mula to the quantum field case, show that it also admits 
a thermodynamic limit and that it satisfies the scaling 
property found for the classical field. We have used the 
quantum approximate formula to evaluate the effect of 
an energy cut-off, not required in the quantum theory 
and unavoidable in the classical field model. 

The perspective of using the condensate phase spread- 
ing to experimentally distinguish among different statis- 
tical ensembles is fascinating, although the measurement 
of the intrinsic phase diffusion of a Bose-Einstein conden- 
sate discussed here remains a challenge and will be the 
subject of further investigations. 
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APPENDIX A: LANDAU-BELIAEV DAMPING 
RATES AND APPROXIMATE PHASE 
DIFFUSION COEFFICIENT IN THE 
THERMODYNAMIC LIMIT 

We start with the Landau damping rate of a Bogoli- 
ubov mode of wave vector q as given in fl3| 



d^kLlj,,{nk - nk')6(eq + ek - ek') (Al) 



with 



Lk,k' = UqVkUk' + {Uq + Vq){UkUk' + VkVk') + VqUkVk'. 

(A2) 

The mode of wave vector q scatters an excitation of wave 
vector k giving rise to an excitation of wave vector k'. 
The final mode has to satisfy momentum conservation so 
that k' = k + q. Energy conservation + ~ ek' is 
ensured by the delta distribution in (jAip . In the integral 
over k we use spherical coordinates of axis q, 9 being the 
polar angle. We introduce the momentum q scaled by 
the inverse of the healing length ^ and the mode energy 
eq scaled by the Gross-Pitaevskii chemical potential pg: 



1/2 



2mpg ^ 

^ = if if 
P9 



2)]i/2. 



(A3) 
(A4) 



As a consequence, the mean occupation number rig is a 
function of q and of the ratio kgT/ pg only, and the mode 
amplitudes Uq, Vq are functions of q only. Introducing the 
notation u = cos 9, one has 



Sieq + ek- ek') = S{u - Ug ) 



2fcg[l + (eg + efe)2]i/2 



where 



l + {eq + ek)'Y^'-{l + e + k') 
2kq 



(A5) 



(A6) 



One can show that Uq is in between —1 and 1 for all val- 
ues of k and q, so that the angular integration is straight- 
forward and leads to 



^k,k' 



q[l + [eq + ekY]"^ 



with 



1-l-fc 



[1 + (eg + efc) 



211/2 



(A7) 



(A8) 



A similar procedure may be applied to the Beliaev 
damping rate for the mode q. From ^is'l one has 



rf = I dPkBlk'{Mk+nk')5{ek + ek'-eq) (A9) 
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with 



with 



Bk,k' = UqUkUk' + {Uq + Vq){VkUk' + UkVk>) + ^^,14 14'. 

(AlO) 

Here the mode of wave vector q decays into an excita- 
tion of wave vector k and an excitation of wave vector 
k'. Momentum conservation imposes k' = q— k. Energy 
conservation e^' = Cq — ek is ensured by the delta dis- 
tribution in (jA9|l . and clearly imposes k < q. With the 
same scaled variables and spherical coordinates as above, 
one obtains 



— (All) 



where 



" 2kq 



(A12) 



One can show that Uq is in between —1 and 1, whatever 
the values of q and k < q, so that angular integration is 
straightforward and gives 



dkBi 



H^q - efe)(l + nk + Uk') 



(A13) 



1 + k 



1 + (eg - 



1/2 



(AM) 



Finally we introduce the rescaled total damping rate, 

(A15) 
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a dimensionless function of ksT/ pg only. From ((48|) one 
then obtains in the thermodynamic limit an approximate 
expression for the phase diffusion coefficient depending 
only on ksT/pg, 



P9 



+ CXD 



dqe {Ai 



(A16) 



with 



_ dk P e^([/fc + Vkfny,{nk + 1) 



< = (f/, + F,)^-e, 



/g °° dkk?- f-lukiuk + 1) 



(A17) 
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